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Abstract: At visible and infrared frequencies, metals show tantalizing 

promise for strong subwavelength resonances, but material loss typically 
dampens the response. We derive fundamental limits to the optical response 
of absorptive systems, bounding the largest enhancements possible given 
intrinsic material losses. Through basic conservation-of-energy principles, 
we derive geometry-independent limits to per-volume absorption and 
scattering rates, and to local-density-of-states enhancements that represent 
the power radiated or expended by a dipole near a material body. We 
provide examples of structures that approach our absorption and scattering 
limits at any frequency; by contrast, we hnd that common “antenna” 
structures fall far short of our radiative LDOS bounds, suggesting the 
possibility for significant further improvement. Underlying the limits is a 
simple metric, IxP/Imj for a material with susceptibility j, that enables 
broad technological evaluation of lossy materials across optical frequencies. 
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1. Introduction 

At optical frequencies, metals present a tradeoff; their conduction electrons enable highly 
subwavelength resonances, but at the expense of potentially significant electron-scattering 
losses [1-10]. In this article we formalize the tradeoff between resonant enhancement and 
loss, deriving limits to the absorption within, the scattering by, and the local density of states 
(LDOS) [11-16] near a lossy, absorptive body of arbitrary shape. Given a material of suscepti¬ 
bility ;|f(cu), the limits depend only on a material enhancement factor \x{0})'\^/lmx{o}) and on 
the incident-beam energy density (leading to a potential 1 /cfi LDOS enhancement for a metal- 



emitter separation d). The power scattered or dissipated by a material body must be smaller 
than the total power it extracts from an incident beam; we show that this statement of energy 
conservation yields limits to the magnitudes of the internal fields and polarization currents that 
control the scattering properties of a body. Unlike previous bounds [6,9,17-27], our limits do 
not depend on shape, size, or topology, nor do they diverge for zero bandwidth. The crucial 
ingredient is that our bounds depend on x are finite only for realistic lossy materials— 
for idealized lossless materials such as perfect conductors, arbitrarily large optical responses 
are possible. We provide examples of structures that approach the theoretical limits, and also 
specific frequency ranges at which common strTictures fall far short. Our bounds apply to any 
absorptive system, and thus provide benchmarks for the response of metals, synthetic plasmonic 
materials (doped semiconductors) [3,5,28-30], and surface-phonon-polariton materials across 
visible and infrared wavelengths, resolving a fundamental question [1-10] about the extent to 
which resonant enhancement can overcome intrinsic dissipation. 

There has been intense interest in exploiting “plasmonic” [2] effects, which arise for mate¬ 
rials with permittivities that have negative real parts, in metals at optical frequencies. Geome¬ 
tries incorporating such materials are capable of supporting highly subwavelength surface res¬ 
onances [1,2]. Yet such a material has inherent loss arising from the typically significant imag¬ 
inary part of %. Even for applications in which absorption is the goal, material loss dampens 
resonant excitations and reduces the overall response. This tradeoff between resonant enhance¬ 
ment and absorption has been investigated for specific geometries amenable to semianalyti- 
cal methods, leading to a variety of geometry-dependent material dependences. For example, 
in the quasistatic limit, coated spheres absorb energy at a rate proportional to Ixl/Imj [31], 
whereas spheroids absorb energy at a rate [32] proportional to Surface modes at pla¬ 

nar metal-insulator interfaces exhibit propagation lengths proportional to (Ree)^ /Ime at very 
low frequencies [33] (in the Sommerfeld-Zenneck regime [2]), but near their surface-plasmon 
frequencies their propagation lengths are approximately proportional to -\/Im£ (cf. App. E). In 
electron energy loss spectroscopy (EELS), the electron scattering cross-section is proportional 
to a “loss function” Im(—1/e) = (Ime)/|ep that enables experimental measurement of bulk 
plasmon frequencies [34,35]. 

We show that |zP/Im;if is a universal criterion for evaluating the optimal response of a 
metal, with a suitable generalization in Sec. 3 for more general media that may be anisotropic, 
magnetic, chiral, or inhomogeneous. For most materials, |z|^/Im;if increases as a function of 
wavelength (as demonstrated in Sec. 3.2), suggesting that if an optimal structure is known, the 
plasmonic response of a metal can potentially be much greater away from its bulk- and surface- 
plasmon frequencies. For effective-medium metamaterials, our bounds apply to both the un¬ 
derlying material parameters as well as to the effective medium parameters, with the smaller 
bound of the two controlling the maximum response. Thus effective-medium approaches can¬ 
not circumvent the bounds arising from their constitutive materials; however, they may find 
practical application if they can achieve resonances at frequencies that are otherwise difficult to 
acheive with the individual materials. 

The limits derived here arise from basic energy considerations. An incident field Einc in¬ 
teracting with a scatterer generates polarization currents P that depend on both the incident 
field and on the shape and susceptibility of the body. A lossy scatterer dissipates energy at a 
rate proportional to the squared magnitude of the currents, |Pp. At the same time, the total 
power extracted from the incident beam, i.e. the “extinction” (absorption plus scattering), is 
proportional to the imaginary part of the overlap integral of the polarization currents with the 
incident field, ^ jy • P, which is known as the electromagnetic optical theorem [32,36-38] 
and can be understood physically as the work done by the incident field to drive the induced 
currents. The overlap integral is only linear in P whereas the absorption depends quadratically 




on P. If the magnitude of P could increase without bound, then, the absorption would become 
greater than extinction, resulting in a physically impossible negative scattered power. Instead, 
there is a limit to the magnitude of the polarization held, and therefore to the scattering prop¬ 
erties of any body comprising the lossy material. We make this argument precise in Sec. 3, 
where we employ variational calculus to derive general limits for a wide class of materials, 
and we also present limits specihc to metals, which are typically homogeneous, isotropic, and 
nonmagnetic at optical frequencies. We consider here only bulk susceptibilities, excluding non¬ 
local or quantum effects [39-41]. The key results are the limits to absorption and scattering 
in Eqs. (29a,29b,31a-32b) and the limits to LDOS enhancement in Eqs. (34a,34b). Before de¬ 
riving the limits, we present the volume-integral expressions for absorption, scattering, and 
radiative and nonradiative LDOS in Sec. 2. 

In Sec. 4 we compare the response of a number of structures towards achieving the various 
limits. Eor far-held absorption and scattering, we hnd that ellipsoidal nanoparticles are ideal 
and can reach the limits across a wide range of frequencies. Eor near-held enhancement of 
power expended by a dipole into radiation or absorption, we hnd that it is much more difficult 
to reach the limits. The nonradiative LDOS near a planar metal surface reaches the limit at 
the “surface-plasmon frequency” of the metal. At lower frequencies, common structures (thin 
hlms, metamaterials) fall far short. At all frequencies, common designs for radiative LDOS 
enhancement exhibit suboptimal response. These results suggest the possibility for signihcant 
design improvement if the limits are achievable (i.e. “tight”). 

Previous limits to the electromagnetic response of metals have emphasized a variety of limit¬ 
ing factors. At RE and millimeter-wave frequencies, where the response can be bounded relative 
to that of a perfect electric conductor (PEC) [42], the Wheeler-Chu-McLean limit [23-25,42] 
bounds the radiative Q factor of an electrically small antenna. At optical frequencies, absorption 
loss increases and often dominates relative to radiative loss. There are known lower bounds on 
the absorptive Q for low-loss, quasistatic structures [26], and more generally for metals of any 
size with susceptibilities comprising Lorentz-Drude oscillator terms [9]. 

Limits to frequency-integrated extinction are also known. Purcell derived the hrst such limit, 
using the Kramers-Kronig relations [38] to bound the integrated response of spheroidal parti¬ 
cles to their electrostatic (o) = 0) induced dipole moments [17]. Recently the limits have been 
extended to arbitrary shapes [18-20], but one obtains a different limit for each shape. More¬ 
over, it is important in many applications to disentangle the single-frequency response and the 
bandwidth, and to do so separately for absorption and scattering. 

Single-frequency absorption and scattering limits have primarily been derived via spherical- 
harmonic decompositions, originally for spherically symmetric scatterers [43,44] and later for 
generic ones [2 1 , 22] . This approach has been generalized recently, yielding limits in terms 
of the inverse of a scattered-held-operator [27], although it the inverse of such an operator is 
seemingly difficult to bound without resorting to spherical harmonics. In Sec. 3 we show that the 
scattered-held-operator approach and our material-dissipation approach share a common origin 
in volume-integral equations. The key distinction is that the scattered-held operator and its 
corresponding limits are independent of material but dependent on structure, whereas our limits 
incorporate material properties and are independent of structure. Both classes of limits apply 
to any linear body. In Sec. 5 we provide a more detailed comparison, hnding that the spherical- 
harmonic limits may provide better design criteria at lower (e.g. rf) frequencies, whereas our 
limits should guide design at higher frequencies, especially in the held of plasmonics. 

This work was partly inspired by our recent bounds [45] on extinction by quasistatic nanos¬ 
tructures. In that work we derived bounds via sum rules of quasistatic surface-integral opera¬ 
tors [46,47]; equivalently, we could have [48] derived the bounds via analogous constraints in 
composite theory [49,50]. The key distinction between this work and our previous work [45] is 
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Fig. 1. Scattering problem under consideration. An incident field Fjne = (Einc ZoHjnc)^ 
impinges on a lossy scatterer with a susceptibility tensor ^(x, to). The shape and topology 
of the scatterer are arbitrary: It may be periodic, extend to infinity, or consist of multiple 
particles. The limits presented in Sec. 3 hinge on the fact that absorption is a quadratic 
functional of the electromagnetic fields, whereas extinction, by the optical theorem, is the 
imaginary part of a linear functional of the fields. In Sec. 3, we present general limits for 
tensor susceptibilities and also simplified limits for metals. 


that here we find limits in the full Maxwell regime, such that our bounds apply to any structure 
at any size scale, and they apply to functions of the scattered fields (e.g. scattered power and 
radiative LDOS), which have zero amplitude in quasistatic electromagnetism. An additional 
benefit of our simplified energy-conservation approach is that we can bound the responses of 
anisotropic, magnetic, and/or inhomogeneous media, whereas the surface-integral sum-rule ap¬ 
proach only works for isotropic and nonmagnetic materials. In this work we also consider the 
local density of states, which we did not consider previously and which represents an important 
design application. 

2. Absorption, scattering, and LDOS expressions 

We consider lossy media interacting with electromagnetic fields incident from fixed external 
sources (e.g. plane waves or dipole sources). Figure 1 illustrates the conceptual setup: a generic 
scatterer with susceptibility tensor x absorbs, scatters, and extinguishes (extinction defined as 
absorbed plus scattered power) incident radiation at rates proportional to volume integrals over 
the scatterer. In this section we present the known volume-integral expressions for absorption 
and scattering, and we also derive volume-integral expressions for the power expended by a 
dipole near such a scatterer, which is either radiated to the far field or absorbed in the near field. 
Relative to a dipole in free space, the enhancement in power expended is given by the relative 
increase in the local density of states (LDOS) [15]. 

The scatterer is taken to consist of a lossy, local material that is possibly inhomogeneous, 
electric, magnetic, anisotropic, or bianostropic (chiral). We assume the scatterer is in vac¬ 
uum, with permittivity Cq, permeability /Tq, impedance Zq = y^/ro/io, and speed of light 
c = 1/ y/eo/fQ. Extending the limits of the following section to non-vacuum and possibily inho¬ 
mogeneous backgrounds is relatively straightforward and is discussed in Sec. 5. The response 
of the scatterer can be described by induced electric and magnetic polarization currents, P(x) 
and M(x), which satisfy the constitutive field relations 

D(x) = eoE(x) -f P(x) 

B(x)=aio[H(x)+M(x)]. 


( 1 ) 







For the general class of materials considered here, the currents P and M each depend on both 
E and H through a unitless 6x6 susceptibility tensor % [15,51]: 

(iivi) ^ (zqh) ^ 

where F is a generalized vector field containing both electric and magnetic fields. For isotropic 
media with relative permittivity Er and relative permeability jir, the susceptibility tensor 
comprises only the diagonal elements — 1 and /Xr — 1. Lossy media as considered here 
have susceptibilities that satisfy the positive-definiteness condition (for an time conven¬ 
tion) [51,52] 

a)lmf>0, (3) 

where Im^ = (^ —X^)/2i (and ^ represents the conjugate transpose). 

When light impinges on the scatterer, the absorbed, scattered, and extinguished powers can be 
written as overlap integrals of the internal currents and fields [53]. The absorption (dissipation) 
within such a medium is the work done by the total fields on the induced currents, given by the 
expression [54]: 


i^abs = ^Im^F*(x).x(x)F(x)dy (4) 

where the asterisk denotes complex conjugation. Equation (4) reduces to the usual 
{eoCo/2) (Imj) fy |Ep for homogeneous, isotropic, nonmagnetic media. 

The total power extracted from the incident fields—the extinction—is the sum of the ab¬ 
sorbed and scattered powers and can be computed by the optical theorem [38,55]. Although 
commonly written as an integral over fictitious effective surface currents [38], the optical the¬ 
orem can also be written as a volume integral over the polarization currents [36,53,56], repre¬ 
senting the work done by the incident fields on the induced currents: 

Fext = ^lmJ^Fl,{x)-^{x)F{x)dV (5) 

where, as for F, we define Fjnc by 

Fine = ) ■ (6) 

V'^oWinc / 

The scattered power is the difference between extinction and absorption: 

Fscat — Fgxt Fabs 

= ^ Im [FL(x) - F*(x)] • f (x)F(x) dV. (7) 

In addition to absorbing and scattering light, structured media can also alter the spontaneous 
emission rates of nearby emitters. Increased spontaneous emission shows exciting potential for 
surface-enhanced Ramam scattering (SERS) [57,58], fluorescent imaging [59,60], thermopho- 
tovoltaics [61,62], and ultrafast light-emitting diodes (LEDs) [63]. The common metric for 
the enhanced emission rate is the (electric) local density of states (LDOS), which represents 
the density of modes weighted by the relative energy density of each mode’s electric field at a 
given position [15]. Equivalently, and more generally, the LDOS enhancement represents the 
enhancement in the total power expended by an electric dipole radiator [15,64,65], into either 


radiation or dissipation. Similar to extinction, the total LDOS is the imaginary part of a linear 
functional of the induced electric helds [12]: 



( 8 ) 


where denotes the held from a dipole source at xq polarized in the S; direction, with a dipole 
moment po = CoS;, and the sum over j accounts for all possible orientations (the conventional 
LDOS corresponds to a randomly oriented dipole [12]). 

To connect the LDOS to the material properties, we rewrite it as a volume integral over the 
helds within the scatterer. The total held at the source position, Ej^. (xq), consists of an incident 
held and a scattered held: 


E.y^(xo) — Einc,5y(xo) T Escat,^^ (xq) ■ 


(9) 


The incident held is known—it is the held of a dipole in vacuum—and can be left as-is (note 
that the imaginary part of a dipole held does not diverge at the source location [66]). The 
scattered held arises from interactions with the scatterer and is the composite held from the 
induced electric and magnetic currents, radiating as if in free space: 
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G^^(xo,x)P(x) + G^''^(xo,x)M(x) 


Escat(xo) — 


( 10 ) 


where and G™ are free-space dyadic Green’s functions [51,67]. We could at this point 


insert Eqs. (9, 10) into Eq. (8) and have a volume-integral equation for the total LDOS. However, 
note that the Green’s functions in Eq. (10) represent the helds of free-space dipoles and the 
excitation in this case is also a dipole. We follow this intuition to replace the Green’s functions 
by the incident helds. 

Whereas the source dipole at xq generates incident helds at points x within the scatterer, the 
Green’s functions in Eq. (10) are the helds at xq from a dipole at x. By reciprocity [68] one can 
switch the source and destination points of vacuum Green’s functions 


Gf (xo,x) = Gf (x,xo) 
Gff{xo,x) = -iioGY{x,xo) 


(lla) 

(llb) 


where for clarity we indexed the Green’s function tensors. Now one can see that the product 
§/ • G^^P equals gQ^'Einc • P. The magnetic Green’s function yields the incident magnetic held, 
with a negative sign arising from reciprocity, as in Eq. (1 lb). Equation (9) can be written 



( 12 ) 


where we have dehned 



(13) 


Inserting the held equation, Eq. (12), into the LDOS equation, Eq. (8), yields the total LDOS 
as the sum of free-space and scattered-held contributions: 



( 14 ) 








where po = aP' is the free-space electric LDOS [69]. It is typically more useful to nor¬ 

malize p to Po; 


plot 

Po 



(15) 


where k — co/c. Equation (15) relates the total LDOS to a volume integral over the scatterer, 
which will enable us to find upper bounds to the response in the next section. 

For many applications it is important to distinguish between the power radiated by the dipole 
into the far-held (where it may be imaged, for example) and the power absorbed in the near 
held (which may productively transfer heat, for example). Absorbed power is given by Eq. (4), 
and thus the nonradiative LDOS enhancement Pnr/po is given by Eq. (4) divided by the power 
radiated by a dipole (of amplitude Co) in free space, fi-ad = EqCO'^/ 12tic^ (Ref. [38]): 



(16) 


Finally, just as the scattered power in Eq. (7) is the difference between extinction and absorp¬ 
tion, the radiative part of the LDOS is the difference between the total and nonradiative parts: 



(17) 


3. Limits 


Given the power and LDOS expressions of the previous section, upper bounds to each quan¬ 
tity can be derived by exploiting the energy conservation ideas discussed in the introduction 
(Sec. 1). The extinction is the imaginary part of a linear function of the polarization currents, 
whereas absorption is proportional to their squared magnitude (and scattered power is the differ¬ 
ence between the two), and thus energy conservation yields finite optimal polarization currents 
and fields for each quantity. 

Just as one can use gradients to find stationary points in finite-dimensional calculus, one 
can use variational derivatives [70] to find stationary points of afunctional (i.e. a function of a 
function). It is sufficient here to consider functionals of the type P = J f*g, which arise in the 
power expressions, Eqs. (4,7,16,17). The variational derivative of P with respect to g is given 
by ^ / f*g = /*, analogous to the gradient in vector calculus: (a^x) = a^. The primary 

distinction is the dimensionality of the space and thus the appropriate choice of inner product. 

The optimal fields for the various response functions P are therefore those for which P is 
stationary under small variations of the field degrees of freedom. The field F is complex-valued, 
such that one could take the real and imaginary parts of F as independent {P is a nonconstant 
real-valued functional and therefore not analytic [71] in F), but a more natural choice is to 
formally treat the field F and its complex-conjugate F* as independent variables [71,72]. Then 
a necessary condition for an extremum of a functional P[F] is for the variational derivatives 
with respect to the field degrees of freedom to equal zero, 5P/5F = 0 and 5P/5F* = 0 (which 
are the Euler-Lagrange equations [70] for functionals that do not depend on the gradients of 
their arguments). Because our response functions are real-valued, the derivatives with respect 
to F and F* are redundant—they are complex-conjugates of each other [71]— and the condition 
for the extremum can be found with the single equation 



(18) 




where we have chosen to vary F* instead of F for its slightly simpler notation going forward. 
We apply this variational calculus approach to bound each response function of interest. First, 
we derive limits for the most general class of materials under consideration. Then we specialize 
to metals, an important class of lossy media that are typically homogeneous, isotropic, and 
nonmagnetic at optical frequencies. 


3.1. General lossy media 

We consider first Eq. (7), for the scattered power. Setting the variational derivative of Pjcat to 0 
yields 

^Fscat £0® 

5¥* Y~ 

The optimal field that satisfies Eq. (19) is 

Fscat,opt(x) = ^ [Imf (x)] f '''(x)Fine(x) (20) 


Z^Finc 
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+ (lmx)F 


= 0 . 


(19) 


for all points x within the scatterer volume V. The optimal field is guaranteed to exist because 
(almx is positive-definite, per Eq. (3), and therefore invertible. We have only shown that this 
is an extremum, not a maximum, but because (olmx is positive-definite, the scattered power in 
Eq. (7) is a concave functional, for which any extremum must be a global maximum [73]. 

One can see that the optimal fields within the scatterer are related to the incident fields (di¬ 
rectly proportional for homogeneous media), which conforms intuitively with the scattered- 
power expression in Eq. (7). The internal field should strongly overlap with the incident field, 
to increase the power extracted from the incident beam, while the susceptibility dependence 
balances between maximizing extinction and minimizing absorption. 

A similar procedure yields the optimal internal fields for maximum absorption within a scat¬ 
terer. Although the absorbed power as given in Eq. (4) is unbounded with respect to F, adding 
the constraint that absorption must be smaller than extinction (i.e. the scattered power must 
be nonnegative) imposes an upper bound. Because Eq. (4) is unbounded, the Karash-Kuhn- 
Tucker (KKT) conditions [74] require that the constraint P^cat > 0 must be active, i.e. P^cat = 0. 
Eollowing standard constrained-optimization theory [74], we define the Lagrange multiplier £ 
and the Lagrangian functional .if = Pabs + ^Fscat- The extrema of .if satisfy 
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( 21 ) 


To simultaneously ensure that the scattered power also equals 0, one can verify that the La¬ 
grange multiplier is given by f = 2. Then the optimal internal fields are 

Fabs,opt(x) = i [imZ(x)] ' Z^(x)Finc(x) (22) 

which are precisely double the optimal scattering fields of Eq. (20). Maximizing Pabs subject 
to Pscat > 0 is a problem of maximizing a convex functional subject to a convex quadratic 
constraint, such that the solution in Eq. (22) must be a global maximum [75]. 

The limits to the scattered and absorbed powers are given by substituting the optimal fields 










in Eq. (20) and Eq. (22) into Eq. (7) and Eq. (4), respectively: 

neat I Kc ■ f' (Imf) ^' f F.ne d\ (23a) 

nbs , nxt I Kc ■ I' (Imf) ^' fF.ne d\ (23b) 

where extinction has the same limit as absorption, which can be derived by maximing Pext 
subject to Pscat > 0. The limits depend only on the intensity of the incident field and the ma¬ 
terial susceptibility x over the volume of the scattered The product x (imz) discussed 

further below, sets the bound on how large the induced currents can be in a dissipative medium. 
Whereas optimal per-volume scattering occurs under a condition of equal absorption and scat¬ 
tering, optimal per-volume absorption occurs in the absence of any scattered power and can be 
larger by a factor of four. This ordering is reversed in the spherical-multipole limits [21,22], 
where the scattering cross-section (not normalized by volume) can be four times larger than the 
absorption cross-section. 

Note that Eqs. (23a,23b) look superficially similar to the absorbed- and scattered-power lim¬ 
its in Ref. [27]. As discussed in Sec. 1, they share a common origin as energy-conservation 
principles applied to integral equations. The key distinction is which quantity serves as a non¬ 
negative quadratic (in F) constraint. We treat absorption as the quadratic quantity, given by 
Eq. (4), with the scattered power as the difference between extinction and absorption. The 
scattered-field-operator approach rewrites the scattered power via a volume integral equation 
(VIE) [51]. This yields a non-negative, quadratic scattered power that is of the same form as 
Eq. (4) except with the replacement Im^ —Im^^, where is a scattered-field integral operator 
with the homogeneous Green’s function as its kernel (the electric component of ^ is used in 
Eq. (10)). Energy conservation leads to the limits of Ref. [27], which take a similar form to 
Eqs. (23a,23b), except Imj —Im^^ and the factors of four are reversed. The limits in Ref. [27] 
are a generalization of the spherical-multipole limits of antenna theory [21,22], which treat the 
special case of the scattered field decomposed into spherical harmonics. 

Our limits have very different characteristics from those of [21,22,27]. Our approach, via 
absorption as the quadratic constraint, yields limits that incorporate the material properties and 
are independent of structure. This naturally results in per-volume limits, a normalization of 
inherent interest to designers. The scattered-field-operator approach yields limits that are inde¬ 
pendent of material but depend on the structure, in a way that can be difficult to be quantified 
because the inverse of the scattered-field operator is not known except for the simplest cases 
(e.g. dipoles). A spherical-harmonic decomposition of the operator yields analytical limits, but 
only to the cross-sections, without normalization. The cross-section is inherently unbounded 
(increasing linearly with the geometric cross-section at large sizes) and thus difficult to use 
from a design perspective. The different normalizations are responsible for the different order¬ 
ings of the absorbed- and scattered-power limits. In Sec. 5 we extend this comparison to show 
that our material-dissipation approach provides better design criteria at optical frequencies. 

The same derivations lead to optimal fields and upper bounds for the radiative and nonradia- 
tive LDOS. The optimal fields are nearly identical in form: 

F,^,rad,opt(x) = ^ (Imf) (24a) 

F.,,nr,opt (x) = / (Imf) ^' f (x) (24b) 

where the complex conjugation arises due to the lack of conjugation in the LDOS expres¬ 
sions (which itself arises because open scattering problems in electromagnetism are complex- 


symmetric rather than Hermitian [15]). Substituting the optimal fields into the LDOS expres¬ 
sions Eqs. (16,17) gives the LDOS limits 
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(25a) 

(25b) 


As for extinction, the limit to the total LDOS is identical to the limit to the nonradiative LDOS, 
which can be proven by maximizing ptot subject to Prad > 0. 

The absorption, scattering, and LDOS limits in Lqs. (23a,23b,25a,25b) depend on the overlap 
integral of the material susceptibility and the incident field. We can simplify the limits further 
by separating the dependencies, which is simple for homogeneous, isoptropic media but can 
also be done for more general media through induced matrix norms [76]. The integrand in 
Lq. (23a) (and each of the other limits) is of the form z^Az, a quantity related to the norm (i.e. 
“magnitude”) of a matrix A. The induced 2-norm of a matrix A, denoted ||A|| 2 , is given by the 
maximum value of the quantity z'Az/z^z for all z 7 ^ 0. The integral in Lq. (23a) can then be 
bounded for general media and arbitrary incident fields: 
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where the dependence on the material susceptibility is now separated from the properties of the 
incident field Fine. The field intensity iFjncP is proportional to the energy density of the incident 
field: 

2^o|Finc| =2 ^o|F™c| 

= Ve ,inc .inc (27) 


where t/g^nc and are the (spatially varying) incident electric and magnetic energy densi¬ 

ties [38]. Generally the incident fields relevant to Pscat and P^bs are beams with nearly constant 
intensity and infinite total energy, for which one should bound the scattered or absorbed power 
per unit volume of material. Given the operator definition and energy-density relation just dis¬ 
cussed, the scattering and absorption limits in Lqs. (23a,23b) simplify: 
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(28b) 


Plane waves are incident fields of general interest. They have equal electric and magnetic 
energy densities and constant intensities /i^c = cUe. me, where c is the speed of light in vacuum. 
The cross-section of a scatterer is defined a = P/Imc representing the effective area the scat- 
terer presents to the plane wave. Because plane waves are constant in space, the absorption and 
scattering bounds are tighter. In the first line of Lq. (26), |Finc|^ can be taken out of the integral, 
which then simplifies to the average value of the norm of j (imz) X- With this modification 



















to Eqs. (28a,28b), the bounds on absorption and scattering cross-sections per unit volume are: 



(29a) 

(29b) 


which apply for general 6x6 electric and magnetic susceptibility tensors. For susceptibilities 
that are only electric or only magnetic, and therefore 3x3 tensors, the bound is smaller by a 
factor of two, since the incident magnetic field cannot drive magnetic currents (or vice versa). 
The LDOS analogue of Eqs. (29a,29b) is not straightfoward, because the incident helds are 
inhomogeneous. Consequently, we leave Eqs. (25a,25b) as the general LDOS limits for inho¬ 
mogeneous media, and derive a simpler version for metals in the next subsection. 

Nanoparticle scattering and absorption are often written in terms of electric/magnetic polariz¬ 
abilities and higher-order moments [22,32,43,77], whereas Eqs. (29a,29b) are bounds in terms 
of only the material susceptibility. One implication is that Eqs. (29a,29b) imply restrictions on 
the number of moments that can be excited, or the strengths of the individual excitations, in a 
lossy scatterer. A lossy scatterer of hnite size cannot have arbitrarily many spherical-multipole 
moments excited, nor can a single scatterer of very small size achieve full coupling to the 
lowest-order electric and magnetic dipole moments. Scatterers for which Imj/jjl ^ 
cannot achieve ^ cross-sections per “channel,” even on resonance. 

3.2. Metals 

Metals represent an important and prevalent example of lossy media. (We define a material to 
behave as a “metal” at a given frequency (O if Rej)®) < — 1, thus including materials such 
as SiC [78,79] and Si02 [79] that support surface-phonon polaritons at infrared wavelengths.) 
At optical frequencies, common metals have homogeneous, isotropic, and nonmagnetic sus¬ 
ceptibilities, enabling us to write the matrix norm of the previous subsection as a simple scalar 
quantity. 
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(30) 


where x{^) the electric susceptibility. Another alteration in the metal case is that the inci¬ 
dent magnetic energy density, Un.mc^ drops out of the limits because the magnetic polarization 
currents are zero in Eq. (2) (and therefore one can simplify Fine,Fine —5" Fine). This is not a 
quasistatic restriction to small objects that only interact with the incident electric field; larger 
objects that potentially interact strongly with the magnetic held remain valid. But their optimal 
response can be written in terms of only the incident electric held, since absorption and extinc¬ 
tion by nonmagnetic objects can also be written only in terms of electric helds, per Eqs. (4,5). 

The limits to per-volume absorption and scattering, simplifying Eqs. (28a,28b), are 
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(31b) 


Similarly, the cross-section limits (reduced by a factor of two relative to Eqs. (29a,29b) because 



















Uh. me is not in the limit) are 
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(32a) 

(32b) 


where as before k= (ojc. Whereas the optimal per-volume scattering occurs at a condition of 
equal absorption and scattering, in App. C we also derive limits under a constraint of suppressed 
absorption, as may be desirable e.g. in a solar cell enhanced by plasmonic scattering [4]. 

The limits to the power expended by a nearby dipole emitter can be similarly simpli¬ 
fied for metals. The incident field is the field of an electric dipole in free space, propor¬ 
tional to the product of the homogeneous Green’s function and the dipole polarization vector. 


Finc.ij = Einc,ij = G^^Sj (because the metal is nonmagnetic, only the incident electric field is 

relevant). The integral over the incident field in Eqs. (25a,25b), summed over dipole orienta- 
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tions, is given by Jv 
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where 


denotes the Frobenius norm [76]. 


For the homogeneous photon Green’s function the squared Frobenius norm is shown in App. G 
to be 
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where the 1 /r^ and \/r^ terms arise from near-field nonradiative evanescent waves, and the 
1 jr^ term corresponds to far-field radiative waves. 

Inserting Fq. (33) into the LOGS limits, Fqs. (25a,25b), yields a complicated integral that 
depends on the exact shape of the body. The integrand is positive, though, so one can instead 
calculate a limit by integrating over a larger space that encloses the body (we will show that 
most of the potential for enhancement occurs very close to the emitter, such that the exact shape 
of the enclosure is usually irrelevant). We consider in detail the case in which the scatterer is 
contained within a half-space, but we also note immediately after Fqs. (34a,34b) the necessary 
coefficient replacement if the enclosure is a spherical shell. All structures separated from an 
emitter must fit into a spherical shell, and thus we have not imposed any structural restrictions 
(in particular, there is no need for a separating plane between the emitter and the scatterer). 

We consider a finite-size approximation of the half-space: a circular cylinder enclosing the 
metal body, a distance d from the emitter and with equal height and radius, L (ultimately we are 
interested in the limit L The volume integrals are straightforward in cylindrical coordi¬ 

nates, yielding j = % j 6d^, fyl/r^ — n/d, and fyl/r^ — n ln(2)L, for d (discarding 
the contributions d/L for the evanescent-wave terms). Then the limits to radiative and non¬ 
radiative LDOS rates are 
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where 0( ) signifies “Big-O” notation [80]. Note that the 0{kL^ terms, which arise from the 
far-field excitation, diverge as the size L of the bounding region goes to oo, whereas one would 


































expect the near-field excitation to be most important. The 0(A:L) divergence as L —>^ oo is un¬ 
physical: it represents a polarization current that is proportional to the 1 /r incident field, ac¬ 
cording to Eqs. (24a,24b), over the entire half-space, maintaining a constant energy flux within 
a lossy medium. Hence, this 0{kL^ term, while a correct upper bound, is overly optimistic, 
and the attainable radiative contribution must be non-diverging in L. One could attempt to sep¬ 
arately bound the evanescent and radiative excitations. However, L also represents the largest 
interaction distances over which polarization currents contribute to the LDOS (in Eq. (10), for 
example); in App. B we show that for reasonable interaction lengths L and near-field separa¬ 
tions d, the contribution of the 0{kL) terms is negligible compared to the \/d^ terms (because 
the divergence is slow). Thus in the near field, where the possibility for LDOS enhancement is 
most significant, the limits are dominated by the 1 /d^ terms: 
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(35a) 

(35b) 


A spherical-shell enclosure of solid angle Q. yields the same result but with the replace¬ 
ment 1/8 Q./47Z in each limit. Again we see the possibility for enhancement proportional 

to |xP/Im;if. There is the additional possibility of near-field enhancement proportional to 
l/{kd)^, which arises from the increased amplitude of the incident field at the metal scatterer. 

Eigure 2 depicts IxP/Imj as a function of wavelength for many natural and synthetic [3, 
5,28-30,79,81-83] metals. Three recent candidates for plasmonic materials in the infrared— 
aluminum-doped ZnO (AZO), and silicon-doped InAs—are included using Drude models of 
recent experimental data from Naik et. al. [29], Law et al. [30], and Sachet et al. [83], respec¬ 
tively. Eor the conventional metals, data from Palik [79] was used; high-quality silver, consis¬ 
tent instead with the data from Johnson and Christy [85] and Wu et al. [86], would have smaller 
losses and a factor of three improvement in IzP/ImZ- ^ broadband version of the metric can 
be computed for extinction or LDOS by evaluation at a single complex frequency [16,56]. 

The material enhancement factor IzP/ItnZ appears in the absorption cross-section of qua¬ 
sistatic ellipsoids [32]; here we have shown that it more generally bounds the scattering re¬ 
sponse of a metal of any shape and size. It arises in the increased amplitude of the induced 
polarization currents; for example, the optimal scattering fields of Eq. (20) simplify in metals 
to the optimal currents 
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with similar expressions for the optimal currents for maximum absorption and LDOS. The feic- 
tor Iz P/ImZ provides a balance between absorption and scattering: in terms of the polarization 
currents, the absorption in a metal is proportional to (Imz/|zP) Iv whereas the extinction 
is proportional to Im -P, thus requiring currents proportional to |zl^/Irnz for absorp¬ 
tion and extinction to have the same order of magnitude. The expression is intuitively appealing 
because a large \x\ signifies the possibility to drive a large current, while a large Imj dissipates 
such a current. Our bounds suggest that epsilon-near-zero materials [87,88], with \x\ « 1, re¬ 
quire a very small Imj to generate scattering or absorption as large as can be achieved with 
more conventional metals. 

A similar, alternative understanding can be attained by considering the currents J = dV/dt = 
—icoP. Defining the complex resistivity of the metal as p = ije^COX, the analogue of Eq. (36) 








Fig. 2. A comparison of the metric \x '\^which limits absorption, scattering, and 
spontaneous emission rate enhancements, for conventional metals (Ag, Al, Au, etc.) [79] as 
well as alternative plasmonic materials including aluminum-doped ZnO (AZO) [29], highly 
doped InAs [30], SiC [79], TiN [81], ITO [82], and Dysprosium-doped cadmium oxide [83] 
(CdO:Dy). Silver, aluminum, and gold are the best materials at visible and near-infrared 
wavelengths, although at higher wavelengths the structural aspect ratios needed to achieve 
the limiting enhancements may not be possible. The dotted lines indicate wavelengths at 
which resonant nanorods would require aspect ratios greater than 30, approximating the 
highest feasible experimental aspect ratios [84]. Despite having lower maximum enhance¬ 
ments, AZO, doped InAs, and SiC should be able to approach optimal enhancements in the 
infrared with realistic aspect ratios. 
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The enhancement factor IzP/Imj thus corresponds to the inverse of the real part of the metal 
resistivity. 


\xi(o)\^ = ^ (38) 
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which corroborates the idea that small metal resistivities enable large field enhancements, as 
discussed recently for circuit [41] and metamaterial [6] models of single-mode response. 

The limits in Eqs. (32a,32b,35a,35b) can be applied additively to multiple bodies; the per- 
volume absorption and scattering limits of Eqs. (32a,32b) are equally valid for a single particle, 
multiple closely spaced particles, layered films, or any other arrangement. Similarly, the LDOS 
limits in Eqs. (35a,35b) can be extended to e.g. a structure confined within two half-spaces, but 
with a prefactor of 2 x 1/8 = 1/4, and any other arrangement in space is similarly possible. 

There are two asymptotic limits in which our bounds diverge: lossless metals (Im^ —0), 
and, for the LDOS bound, the limit as the emitter-metal separation distance (7 —> 0. In each case, 
the divergence is required, as there are structures that exhibit arbitrarily large responses. Eor 






example, as the loss rate of a small metal particle goes to zero, it is known that the absorption 
per unit volume increases until, for a given size, the radiative loss rate equals the absorptive 
loss rate [43,44,89,90]. However, if the size (and therefore the radiation loss rate) is decreased 
concurrently with the material loss, the cross-section per unit volume can be made arbitrarily 
large. Thus, for a small enough particle, any (Jext/T is possible, and the limit must diverge 
as material loss approaches zero (regularized physically by both nonzero loss and nonlocal 
polarization effects [39^1]). Similarly, the LDOS can diverge in the limit of zero emitter- 
scatterer separation, both for lossy materials where absorption diverges [91] and for lossless 
materials with sharp corners [92]. The latter case can be reasoned as follows: the helds at a sharp 
tip, either dielectric or metal, diverge for any nonzero source (of compatible polarization) [38, 
93,94]. By reciprocity, for a source inhnitesimally close to the tip, the LDOS must diverge. 

4. Optimal and non-optimal structures 

We turn now to the design problem: are there structures that approach the limiting responses 
set forth by Eqs. (32a,32b,35a,35b)? We show that optimal ellipsoids can approach both the 
absorption and scattering limits across many frequencies by tuning their aspect ratios. For the 
LDOS limit, however, the optimal designs are not as clear. At the resonant (“surface-plasmon”) 
frequency Osp of a given material the prototypical planar surface exhibits a nonradiative LDOS 
enhancement approaching the limit of Eq. (35b). However, at lower frequencies, neither thin 
hlms [95] nor common metamaterial approaches for tuning the resonant frequency achieve the 
|Xp/Im;if enhancement, thereby falling short of the limit. Similarly, representative designs for 
increased radiative LDOS are shown to fall orders of magnitude short of the limits. These struc¬ 
tures fall short because the near-held source excites higher-order, non-optimal “dark” modes 
that reduce the LDOS enhancement. 

To compute the electromagnetic response of the structures in Fig. 3-5, we employed a free- 
software implementation [96,97] of the boundary element method (BEM) [98]. Where possible 
(quasistatic ellipsoid extinction [32], planar metal LDOS [99]), we also used exact analytical 
and semi-analytical results. 

4.1. Absorption and scattering 

Small ellipsoids approximated by their dipolar response can approach the limits of Sec. 3 across 
a wide frequency range by tuning their aspect ratios. Ellipsoids reach the absorption limits 
for small (ideally quasistatic) structures, and the scattering limits for larger structures that are 
still dominated by their electric dipole moment. The quasistatic absorption cross-section of an 
ellipsoid, for a plane wave polarized along one of the ellipsoid’s axes, is [32] 



(39) 


where L is the “depolarization factor” (a complicated function of the aspect ratio) [32] along 
the axis of plane wave polarization. The optimal response is achieved for the aspect ratio such 
that L = Re [—\/x{(i>)), which yields a polarization held within the particle of [32]: 



(40) 


which is exactly the optimal absorption condition of Eq. (22), as can be seen by comparison 
with Eq. (36). For this optimal depolarization factor, the peak absorption cross-section per unit 





Fig. 3. Absorption (blue) and scattering (red) cross-sections per unit particle volume for 
nanoparticles of (a) gold [79] and (b) Si-doped InAs [30], illuminated by plane waves 
polarized along the particle rotation axis. The ellipsoid aspect ratios can be tuned to ap¬ 
proach both the maximum absorption and maximum scattering cross-sections (black) of 
Eqs. (32a,32b). The dimensions of the nanoparticles are optimized at three representative 
wavelengths, in the visible for gold (constrained to have radii not less than 5nm) and at 
longer infrared wavelengths for doped InAs. Whereas the maximum-ahsorption particles 
are small to exhibit quasistatic behavior, represented by dashed lines in (a), the maximum¬ 
scattering particles are larger such that their scattering and absorption rates are equal. 


volume is given by 

f7abs(dl) 

V 

thereby reaching the general limit given by Eq. (32b). Equation (41) is valid for both oblate 
(disk) and prolate (needle) ellipsoids. Here we have considered the cross-section for a single 
incident plane wave; if one were interested in averaging the cross-section over all plane-wave 
angles and polarizations (as appropriate for randomly oriented particles), then it is possible to 
find a bound that is tighter, by 33% for most materials, in the quasistatic regime. In that case 
the bounds are achieved by disks but not needles. 

Whereas the absorption cross-section is maximized for very small particles approaching the 
quasistatic limit—necessary to exhibit zero scattered power, a prerequisite for reaching the ab¬ 
sorption bounds—the optimal scattering cross-section is achieved for larger, non-quasistatic 
particles that couple equally to radiation and absorption channels. One can show either through 
a modified long-wavelength approximation [100-102] or by coupled-mode theory [43,44] that 
the dimensions of a small particle can be tuned such that the absorption and scattering cross- 
sections are equal, at which point the scattering cross-section per volume is a factor of four 
smaller than the optimal quasistatic absorption. We validate this result with computational op¬ 
timization and show that it enables the design of metallic nanorods with nearly optimal perfor¬ 
mance. 

Figure 3 shows the per-volume absorption and scattering cross-sections of (a) gold and 
(b) Si-doped InAs [30] nanorods designed for maximum response across tunable frequen¬ 
cies. As in Ref. [30] we employed a Drude model for the doped InAs, with plasma frequency 
COp = 2nc/5.5pLm and damping coefficient 7 « O.OSScUp, as is appropriate for a doping density 
on the order of lO^^’cm^^ (Ref. [103]). We employed a free-software implementation [104] 
of the controlled random search [105,106] optimization algorithm to find globally optimal el- 


\x{(o)r 


— _ 

ellinsnid Im;7(a))’ 


(41) 
























Metal 


Al Ag Au AZO InAsSiC 



-LDOS limit, Eq. (20) 

-LDOS, computed 



Emitter-Metal Separation Distance, d 


Fig. 4. Nonradiative LDOS enhancement for randomly oriented dipoles above a flat bulk 
metal, (a) Enhancement as a function of wavelength. For each metal except gold—which 
has significant losses—the nonradiative LDOS at the surface-plasmon frequency ft)sp (dot¬ 
ted line) approaches the limit given by Eq. (35b). The emitter-metal separation distance is 
fixed at d = O.lc/cOsp- The limits are equally attainable for conventional metals such as Al 
and Ag as for synthetic metals such as AZO [29] and highly doped InAs [30], and for SiC. 
(b) Enhancement as a function of metal-emitter separation distance d, with the frequency 
fixed at the surface-plasmon frequency cUsp for each metal. The limiting enhancements are 
asymptotically approached as the separation distance is decreased, because the quasistatic 
approximation of Eq. (42) becomes increasingly accurate. 


lipsoid radii. For the gold nanorods a minimum radius of 5nm was imposed as representative 
of experimental feasibility [39,107] and a size scale at which nonlocal effects are expected 
to remain small [39^1,108]). The gold particles optimized for absorption fall slightly short 
of the limits due to the minimum-radius constraint; in the quasistatic limit (dashed), the ab¬ 
sorption cross-section reaches the limit, as expected from Eq. (41). Both gold and doped-InAs 
nanoparticles closely approach the scattering and absoption limits of Eqs. (32a,32b). Eor Drude 
models, the factor jYmx that appears in both limits simplifies to (sy^jjc, a material con¬ 
stant independent of wavelength, clearly seen in Eig. 3(b). The increase in ff/V as a function 
of wavelength in Eig. 3(a), for gold, can be seen as a measure of the deviation of the material 
response [79] from a Drude model. A constant response for Drude models is not universal; 
near-field interactions, in particular the local density of states (LDOS), instead depend only 
on IzP/Imz ~ CUp/yo), thereby increasing at longer wavelengths, away from the bulk and 
flat-surface plasmon frequencies. 

Not all small particles reach the limiting cross-sections. Coated spheres are common struc¬ 
tures for photothermal applications [31, 109, 110], but their absorption cross-section per unit 
particle volume is proportional to (2/3)|z|/ImZ instead of |zP/Imz (cf. App. D). Their en¬ 
hancement does not scale proportional to \x\^ due to the small metal volume fractions required 
to tune the resonant frequency. 

4.2. LDOS 

Designing optimal structures for the local density of states (LDOS) enhancement limits, 
Eqs. (34b,34a), is not as straightforward as designing optimal particles for plane-wave absorp- 


























tion or scattering. Because the structure is typically in the near field of the emitter, it is difficult 
to design a resonant mode that exactly matches the rapidly varying field profile of the emit¬ 
ter. We show that it is possible to reach the nonradiative LDOS limits at the surface-plasmon 
frequency of a given metal, but that away from these frequencies typical structures fall short. 
Similarly, for the radiative LDOS limits, common structures fall short of the limits, especially 
at longer wavelengths. 

Planar layered structures support bound surface plasmons that do not couple to radiation, and 
thus only improve the nonradiative LDOS. As discussed in Sec. 1, this is potentially useful for 
radiative heat transfer applications, where near-held emission and absorption have been exten¬ 
sively studied [111-115]. Morover, adding either periodic gratings or even random textures can 
couple the bound modes to the far held [116-118] and potentially result in substantial increases 
to the radiative LDOS. Thus we hrst study how closely surface modes in planar structures can 
approach the nonradiative limits, and then we analyze the performance of representative cone- 
and cylindrical-antenna structures relative to the radiative LDOS limits. 

At the surface-plasmon frequency, the prototypical metal-semiconductor interface that sup¬ 
ports a surface plasmon exhibits a nonradiative LDOS approaching the limiting value of 
Eq. (34b). In the small-separation limit (kd <C 1) and at the surface-plasmon frequency cUsp, 
the local density of states near a planar metal interface reduces to [12] 

Pnr (t^sp) 

Po 

thereby approaching exactly the nonradiative LDOS limit of Eq. (34b). (Note that we de- 
hne Po as the electric-only free-space LDOS, different by a factor of two from the elec- 
tricH-magnetic LDOS in [12].) Although the surface-plasmon frequency is typically defined [2] 
as the frequency at which Ree = —1, this is the frequency of optimal response only in 
the zero-loss limit. More generally, we define the surface-plasmon frequency Osp such that 
Re^ (Osp) = Re (—l/x(a)sp)) = 1 /2. Eor gold, which never satisfies Re^ = 1/2 due to its high 
losses, we define surface-plasmon wavelength to be Ajp = 510nm, where Re^(a)) is a maxi¬ 
mum. 

Eigure 4(a) compares semianalytical computations of the nonradiative LDOS near a flat, 
planar metallic interface to the nonradiative LDOS limits given by Eq. (34b). Six metals are 
included: A1 (black), Ag (blue), Au (red), AZO (green), InAs (teal), and SiC (purple), with 
the surface-plasmon frequency Wsp of each in a dotted line and a fixed emitter-metal spacing 
of d = O.lc/cUsp. Every metal except gold—which is too lossy—reaches its respective limit; 
it is possible that a different, nonplanar gold structure, with the correct “depolarization fac¬ 
tor” (VIE eigenvalue, cf. App. A), could approach the limit. Eigure 4(b) shows the emitter- 
metal separation-distance dependence for O) = cUsp. The limits are approached—again, except 
for gold—as the emitter-metal separation decreases and the approximate 1 /d'^ dependence of 
Eq. (42) becomes more accurate. 

There are a few common approaches to tune the resonant frequency below Osp. A standard 
approach is to use a thin film [2, 119], coupling the front- and rear-surface plasmons to cre¬ 
ate lower- and higher-frequency resonances. Other approaches include highly subwavelength 
structuring, to create hyperbolic [120-122] or elliptical metamaterials with reduced effective 
susceptibilies. We show here that such structures do not exhibit the material enhancement fac- 
tor, |zP/IniZ> thus do not approach the limit to Pm given by Eq. (34b). 

The nonradiative LDOS of a thin film can be computed by decomposing the dipole excitation 
into plane waves (including evanescent waves), which reflect from the layers according to the 
usual Eresnel coefficients. The LDOS near a thin film is well-known as an integral over the 
surface-parallel wavevector [69]; for a dipole with fixed frequency (O and height d above the 
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Fig. 5. Away from the surface-plasmon frequency of a given metal—taken here to be 
silver [79] —it is more difficult to reach the radiative and nonradiative LDOS limits, 
Eqs. (34b,34a). The emitter-metal separation d is fixed at li = lOnm for (a) and (b). (a) 
Nonradiative LDOS enhancements, Pnr/Po, for thin films (red) for various silver thick¬ 
nesses (r), and for (type-I) hyperbolic metamaterials (FIMMs, purple) for two silver fill 
fractions (ff). (b) Radiative LDOS enhancements, Pi-ad/po for cone and cylinder antennas, 
with dimensions optimized at wavelengths from A = 450nm to A = 850nm. (c) Scaling of 
Pm- and Pmd for optimized thin films and cone antennas, respectively, as a function of d 
(inset: log-log scale). The scaling of the optimal design appears to be l/d?, with the struc¬ 
tures falling short of their respective limits [the dashed line is (Pi-ad/Po)max] because they 
do not exhibit a |;|;|^/Im;f enhancement. 


film, and a film of optimal thickness f, one can show that Pm is (cf. App. F) 

' Pnr(ftl) 

. Po 



(43) 


which is valid for O) << Wsp (otherwise a bulk is optimal and Eq. (42) describes the response) 
and for relatively small loss, Imx <C |Rex| (as is typical at optical frequencies), to ensure the 
large-wavevector modes are resolved. Unlike a planar interface, an optimal thin film does not 
exhibit the lAp/ImA material enhancement. The thin film falls short because it relies on near- 
held interference to couple the front- and rear-surface plasmons, yielding a resonance that cou¬ 
ples strongly to the dipole emitter over only a small bandwidth of wavevectors (Akp Im;if / |;k |) 
that cancels the resonant enhancement from decreased loss. An alternative understanding arises 
from viewing a thin film as the single-unit-cell limit of a layered hyperbolic metamaterial 
(HMM) [95]. Hyperbolic metamaterials exhibit anisotropic effective susceptibilities such that 
the resonant frequency can be tuned, but their resonances occur within the bulk rather than 
along the surface, such that they cannot yield infinite LDOS even in the limit of zero loss. In 
Ref. [95] we show that the LDOS near an optimal HMM is nearly identical to the LDOS near 
an optimal thin film, as verified in Fig. 5. 

Achieving the radiative LDOS limit is a similarly challenging problem. To reach the radia¬ 
tive LDOS limits, the polarization field must exhibit the Xjr^ spatial dependence of the incident 
field, while also coupling to far-field radiation channels. Here we consider two representative 
structures for tuning the LDOS resonant frequency: mirror-image cones (akin to bowtie anten¬ 
nas) and cylindrical antennas (scaled [66] shorter than A/2), and we show that each falls short 
of the limits across optical frequencies. 

Figure 5 shows the nonradiative and radiative LDOS near silver thin films, effective-medium 
hyperbolic metamaterials (HMMs), cones, and cylinders. In Fig. 5(a) and Fig. 5(b), the emitter- 
metal separation distance is fixed at c/ = lOnm, and the structures are optimized at four wave¬ 
lengths, A = [450,600,725,850]nm, using a standard local optimization algorithm [123]. The 
optimal lengths and thicknesses are included in the figure. The optimal cone half-angles and 
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Fig. 6. A schematic comparison of absorption and scattering limits: multipole limits [21, 
22] to the total cross-section can provide design guidelines at low frequencies, where it 
is difficult to achieve “plasmonic” resonances, but at higher frequencies our dissipation- 
based limits provide tighter limits to the per-volume response. The frequencies at which 
our bounds can be reached range from the bulk plasma frequency (ft) = cOp) down to ft) ~ 
(Op/ARmsx, where A/?max is the maximum achievable aspect ratio. Included is the relevant 
range for silver ellipsoids (red text), assuming = 30. Plasmonic behavior at longer 

wavelengths is possible with materials such as AZO and doped InAs (cf. Fig. 2). 


cylinder radii vary slightly but are typically on the order of 0 = 20°and r = lOnm, resp. Away 
from Osp, the optimal nonradiative LDOS of the thin hlms in Fig. 5(a) is very well described 
by Eq. (43), which shows that the structures fall short of the limits by the material enhancement 
factor. A similar effect is seen for the radiative LDOS in Fig. 5(b), as the cones also fall increas¬ 
ingly short of the limits as the frequency is decreased (« 250 x at A = 850nm). The quantum 
efficiencies (Prad/Ptot) of the cone antenna varies from 50-70%, near the optimal ratio of 50% 
to reach our bounds, suggesting that the reason they fall short is due to a mismatch between the 
emitter and resonance field profiles, not due to coupling to nonradiative channels. Figure 5(c) 
shows the LDOS dependencies as a function of c/ for A = 600nm, in linear (main) and loga¬ 
rithmic (inset) scale, with the structure parameters optimized for each d. The optimal structures 
appear to exhibit the 1/d^ scaling of the limits in Eqs. (34a,34b). An important open question 
is the extent to which structures can be designed to approach the limits, thereby improving over 
current designs by two to three orders of magnitude. 

5. Extensions and discussion 

The limits derived in Secs. (2,3) apply to a general class of materials embedded in vacuum, 
without any scatterers. They can be generalized for non-vacuum backgrounds. For a metal in 
a homogeneous, lossless background with permittivity gbg, the limits are of the same form but 
with replacements k -A ^/£^k and j —?> (£(a)) — ebg)/eba. A similar generalization applies for 
general lossy media in non-vacuum backgrounds. If there are background scatterers present 
(possibly periodic [124-126]), then the derivation is identical except that the Green’s func¬ 
tions are the Green’s functions in the presence of the background scatterers, and the “incident” 
helds therefore incorporate the effects of the background scatterers. The minimum thickness 
of a metal absorber on a substrate [127,128], for example, could be computed with Eq. (32b), 
with the replacement Ejnc = Eq (for substrate reflectivity r), or more generally by 

bounding the incident held via jisinc I < 2 |£’o|- 













Previous approaches to general limits via energy conservation have bounded the response of a 
structure via its scattered-field operator [27]. This has yielded limits to the cross-section, a, for 
spherically symmetric [43,44] or more general [21,22] scatterers whose response has been de¬ 
composed into spherical multipoles. The cross-section limits are proportional to {N^ +A^), 
where N is the number of excited multipoles. 

From a design perspective, there are serious impediments to using such cross-section limits 
(as opposed to our a/V limits). The cross-section itself is unbounded, increasing with the 
size of a large particle [32]. Thus implicit in any use of such bounds is a size normalization, 
but this can only be done at very small scales, where the number of multipoles is 1 (or 2 
for perfect conductors) and our bounds are tighter, and at very large size scales, approaching 
the geometric-optics regime. At intermediate sizes, it is difficult to estimate the number of 
multipoles without further modeling of a given structure. We have presented a new approach, 
using material dissipation, instead of the scattering operator, as the binding constraint. Our 
limits have the unique feature that they incorporate the material properties and are independent 
of structure. From a design viewpoint this is a significant advantage, since for any problem 
there are infinitely many possible structures but typically only a handful of relevent materials. 
Furthermore, the normalization to geometric volume, e.g. G/V, emerges naturally in our limits. 

Reaching our absorption and scattering limits likely requires significant polarization currents 
throughout the scatterer volume, as can be seen in the optimal field profiles in Sec. 3. At low 
frequencies, it is difficult to fabricate structures with sizes or aspect ratios necessary to achieve 
such resonances. For a Drude-metal (appropriate at low frequencies) ellipsoidal nanorod, the 
optimal aspect ratio [32] for maximum absorption scales as (ojcOp, where (Op is the bulk plasma 
frequency. Thus, if we define a maximum feasible aspect ratio ARmax, the minimum frequency 
at which a plasmonic resonance can be achieved is proportional to Op/ARmax- Below this fre¬ 
quency, the multipole limits can serve as a design guide, although the uncertainty about the 
number of multipoles, and the potential mismatch between a non-spherical object and its bound¬ 
ing sphere, remain as barriers. Figure 6 depicts schematically depicts which bound provides 
better design criteria as a function of frequency. Included in Fig. 6 is red text corresponding 
to the limiting frequencies at which silver nanorods (using experimental material data [79] and 
making no Drude approximation) can approach our limits, assuming a realistic [84] maximum 
aspect ratio of 30. 

Similarly, it may be difficult to reach our limits with larger, wavelength-scale solid particles 
that are much larger than the skin depth. One of the conclusions from our work is that such 
particles are particularly inefficient scatterers, and thus should be avoided, because currents 
cannot be excited throughout such a large portion of their volume. At optical frequencies, any 
technology must ultimately incorporate some collection of ordered or disordered scatterers, 
whether in planar arrays [4,57, 129-132], aqueous environments [133-135], or some other 
configuration. Thus even if the individual scatterers have small cross-sections, there can be 
many of them (due to their small volumes), providing a large collective cross-section [107] 
while maintaining the per-volume response of the individual scatterers. 

Aside from particle scattering, our limits extend to situations that do not have multipole 
counterparts. They yield meaningful limits for extended structures (whose large size would 
excite many spherical harmonics), and for the LDOS (where a near-held source would excite 
many spherical harmonics). 

An interesting aspect of our limits is the fact that they apply to any open scattering problem. 
Open systems—in which energy can enter and exit—are typically described by non-Hermitian 
operators. Non-Hermitian operators are not guaranteed to be diagonalizable, and therefore may 
not have a complete basis of eigenfunctions (technically, even for Hermitian operators, rigor¬ 
ous eigendecomposition in inhnite-dimensional spaces is subtle and subject to obscure counter- 


examples [136,137]). Breakdown of diagonalizability only occurs at “exceptional points” that 
must be forced [138,139] and which occur by chance with zero probability. Near an excep¬ 
tional point, however, it is possible to have eigenfunctions that are nearly “self-orthogonal,” 
with exceptionally large modal overlap (or in theory highly nonnormal and ill-conditioned 
Maxwell operators [140]) leading to effects such as destructive interference in scattering “dark 
states” [141] (e.g. Fano resonances [142]) and the Petermann factor for noise enhancement in 
lasers [143,144]. Nevertheless, passive systems near such exceptional points cannot exceed our 
limits, which impose only conservation of energy. 

The limits presented here suggest new design opportunities with metals. Nearly lossless met¬ 
als [8] could manifest unprecedented responses. Even for conventional lossy metals, large-area 
structures that achieved the absorption or scattering limits presented here could potentially do 
so with thicknesses approaching a single atomic layer. Nonlocal interactions, for which it may 
not be possible to separate material and structural properties, would be important in such a 
structure. Finding limits incorporating nonlocal effects would represent an important extension 
to this work. Designing structures to approach the LDOS limits could impact applications such 
as imaging, where there are potentially orders of magnitude improvement to be gained. We 
derive limits for the problem of near-held radiative heat transfer, where the sources are embed¬ 
ded within the designable media, in an upcoming publication [145], and it would be interesting 
to extend the limits derived here to other hgures of merit, potentially hnding new metrics or 
structures for optimal light-matter interactions. 

A. Alternative understanding of the limits: VIE approach 

Here we present an alternative viewpoint for understanding our limits, as arising from sum 
rules over eigenmodes of the volume integral equations of electromagnetism. This approach 
appears to only work for materials with a scalar x (either electric or magnetic), and thus is 
less general than the derivation in the text (we assume a nonmagnetic medium). We include 
this appendix because higher-order sum rules may yield tighter limits in certain scenarios, e.g. 
angle-averaged incident helds. First we show how limits arise from eigenmodes of the volume 
integral equations (VIEs), which can be considered “material resonances.” This connection 
was partially recognized by Rahola [146], and may be related to eigenvalues in SALT laser 
theory [147], but has not since been pursued any further. Material resonances are common in 
quasistatic electromagnetism [148,149], where there is no frequency; here, we show how to 
extend the concept to hxed, nonzero frequencies. 

Inherent to the concept of a resonance in physics is the resonant frequency: intuitively, the 
frequency at which an electromagnetic, elastic, quantum mechanical, or any other type of wave 
oscillates without external forcing in a specihc, predehned structure. In a closed or periodic 
structure, these resonances correspond mathematically to eigenvalues of the underlying differ¬ 
ential equations. Eor photonic structures, defined by a spatially dependent permittivity e(x), 
resonant frequencies o,, are eigenvalues of the eigenequation 



(44) 


dehned with appropriate boundary conditions and the divergence equation V • eE = 0. For a 
material with negligible dispersion (metals introduce further complications), the Maxwell op¬ 
erator ^ = e^*V xVx is frequency-independent, such that the eigenvalue corresponds to a 
resonant frequency. In open systems, resonances are complex poles of the Green’s function (or 
of the scattering operator) rather than true eigenvalues, because their “eigenfunctions” (“leaky 
modes”) diverge exponentially in space [150,151]. For a sign convention of time depen¬ 
dence, these poles (0„ must lie in the lower-half of the complex-frequency plane as shown 
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Fig. 7. A comparison of the (a) resonant-frequency and (b) resonant-susceptibility frame¬ 
works in electromagnetism. The conventional resonant-frequency approach is depicted in 

(a) : the operating frequency is real-valued and can in theory be approached arbitrarily 
closely by a resonance with small imaginary part —lmco„ (i.e. a high-Q resonance). Con¬ 
versely, volume integral equations yield the resonant-susceptibility approach depicted in 

(b) : metal losses, which correspond to Imf = Im(—1/;|^) > 0, inherently impose a mini¬ 
mum separation q to how closely a material resonance—restricted to lie on or below the 
real line—can approach the real system parameters. Moreover, quasistatic structures have 
real-valued eigenvalues, and thus have the potential to achieve the minimum eigenvalue 
separation and maximum optical response. 

in Fig. 7(a). Frequency resonances in electromagnetism are well understood, but material 
resonances—which arise as eigenvalues in integral equations—have hardly been explored at 


all. 


The electric field integral equation (EFIE) formulation of Maxwell’s equations is derived 
through the use of Green’s functions [51,152]. As depicted in Eig. 1, we consider generic scat¬ 
tering problems in which a structure with susceptibility % interacts with an externally imposed 
incident field Ei„c. The response of the scatterer is given by a convolution of the free-space 
Green’s function G with the induced polarization currents P = x^- As in the text we assume 
a scatterer embedded in vacuum, with straightforward generalizations. The total field E is the 
sum of the incident and scattered fields [51,153], 



(45) 


for all points in space, where we choose a negative sign convention for the Green’s function 
(opposite that of Eq. (10) in the text). Equation (45) can be desingularized [154], but our treat¬ 
ment depends only on an abstract spectral decomposition that does not require us to grapple 
with such details. A similar integral equation arises in quantum mechanics, where it is known 
as the Lippmann-Schwinger equation and the susceptibility is replaced by the scattering poten¬ 
tial [155]. Eor a scatterer with homogeneous susceptibility, x is constant and can be taken out 
of the integrand in Eq. (45). 

Homogeneous scatterers are thus defined by the single material parameter X- Resonances of 
Eq. (45), even in open systems, are true eigenvalues because the integral equation has unknowns 
E(x') defined over the finite scatterer domain V, and thus the corresponding eigenfunctions are 
normalizable. Eigenfunctions E„ and eigenvalues of the Green’s function integral operator 
satisfy: 
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for all points x in V. Given the eigenvalue if we choose a material X — Xn = — ^then by 
comparison with Eq. (45) we see that a “standing-wave” E 7 ^ 0 is possible even for Einc = 0. 
Equation (46) is the integral-equation analogue of Eq. (44) (for a homogeneous scatterer), and 
yet we see that the integral operator on the left-hand side of Eq. (46) depends on the structure 
and the frequency but not on the susceptibility. Instead, the eigenvalue of the mode yields 
a resonant value Xn —a resonant material susceptibility, for a fixed frequency. Just as “leaky 
modes” in Eq. (44) are not actually physical solutions of Maxwell’s equations, the resonances 
Xn are not physically realizable materials; we will see below that ImXn < 0 for O) > 0 , as shown 
in Eig. 7(b), corresponding to gain required to overcome modal radiation loss. 

Eor both frequency and material resonances, the separation between the system parame¬ 
ter (e.g. operational frequency) and the resonance defines the magnitude of the response— 
the smaller the separation, the larger the response. In the resonant-frequency framework it 
is difficult to provide a lower bound on the imaginary part of the resonant frequency (0„ 
(thereby bounding the maximum Q). In the resonant-material approach, however, the system 
parameter—the susceptibility, instead of the frequency—has a nonzero imaginary part for a 
lossy system. By causality [156] (or passivity [52]), frequency resonances for a fixed structure 
and material lie in the lower half of the complex-O) plane. Similarly, for a fixed frequency O) > 0, 
the material resonances must reside in the lower half of the complex-^ plane (otherwise, one 
could construct a passive linear material that violates the condition on the resonant frequen¬ 
cies [157]). Thus, as depicted in Eig. 7, there is a minimum separation q = Im^(a)) between 
the material parameter ^(w) = — l/x(a)) and the eigenvalue There are further benefits to 
the resonant-material approach: the solutions to the integral equation are defined only over the 
scatterers, rather than all space, and quantities like the extinguished power and the local density 
of states can be written as volume integrals, ideally suited to a VIE framework. 

To simplify further analysis, we rewrite the VIE of Eq. (45) as: 

(,^ + X^)e = einc (47) 

where the fields are e and Cinc (with lower-case e denoting vector fields restricted to the vol¬ 
ume y, forming a Hilbert space, as opposed to the fields E defined everywhere in space), 
is the identity operator, and ^ is the Green’s function integral operator defined by = 

G(x,x', a))E(x'). The eigenfunctions of ^ are solutions of Eq. (46) and are given in vector 
notation by 

= (48) 

where n is the mode index, = — IjXn, Imi*„ < 0. There are advantages to considering 
the integral operator ^ rather than the Maxwell operator ^ =(l/e)VxVx. Eor the operator 
metals are difficult to treat: material dispersion renders the eigenvalue problem nonlinear 
in 0 )^, and material loss yields a non-Hermitian operator even for closed or periodic struc¬ 
tures [158-160]. The Green’s function operator avoids these difficulties because it assumes a 
fixed frequency and is independent of material, and therefore of material loss. An additional 
advantage of the VIE approach is the compact domain, which sidesteps the subtle normal¬ 
ization [159,161-164] required in the resonant-frequency approach, where the “leaky” fields 
diverge as they extend to infinity. 

We consider open systems—in which energy can enter and exit—such that typical operators, 
including the curl-curl operator ^ of Eq. (44) and the Green’s function integral operator 
of Eq. (47), are not Hermitian. By reciprocity is complex-symmetric [11], such that if it is 
diagonalizable, its generic eigendecomposition can be written [165] 


(49) 


where E is a diagonal operator with entries and is the basis of eigenfunctions e„. 
The assumption of diagonalizability (the existence of a “spectral” eigendecomposition of 
the operator) is commonplace is physics. Technically, even for Hermitian operators, rigorous 
eigendecomposition in infinite-dimensional spaces is subtle and subject to obscure counter¬ 
examples [136, 137]. However, if one makes the reasonable conjecture that the system can 
be simulated on a computer, i.e. that one has a convergent finite-dimensional discretization, 
then one should be able to apply the eigendecomposition in this finite-dimensional system to 
arbitrary accuracy; this can be viewed as a justification of the commonplace assumption of 
diagonalizability. Even in a finite-dimensional problem, of course, diagonalizability of a non- 
Hermitian matrix is not guaranteed [165], but breakdown of diagonalizability only occurs at 
“exceptional points” that must be forced [138,139]. Exceptional points occur by chance with 
zero probability; near an exceptional point, the operator is highly nonnormal [140] and ill- 
conditioned (i.e. cond('^'^^) ^ 1 ), but it is diagonalizable, and we show in the text that the 
response of such structures cannot surpass our limits. By continuity our limits also apply at 
non-diagonalizable exceptional points, thereby justifying our assumption of diagonalizabil¬ 
ity. Our decomposition is similar in spirit to the singularity and eigenmode expansion meth¬ 
ods, SEM and EEM, respectively, for surface integral equations in electromagnetic scattering 
[68,166-169]. 

The complex-symmetry of the operator leads to an atypical (unconjugated, indefinite) “in¬ 
ner product” under which the eigenfunctions are orthogonal; efey = /yE, Ej = 5,;. Because 
modes are not orthogonal under the typical inner product e^e = /y E* • E, they are not power- 
orthogonal [162]: it is possible for energy in one mode to mix into another, leading to effects 
such as destructive interference in scattering “dark states” [141] (e.g. Eano resonances [142]) 
and the Petermann factor for noise enhancement in lasers [143,144]. Another possibility is that 
fyEi E, = 0, i.e. “self-orthogonality,” which renders the eigenfunction basis incomplete [170] 
and corresponds to exceptional points [139, 171]. As shown in the text, energy conservation 
prevents any of these phenomena from surpassing our limits. 

A modal decomposition of the electric field e follows from the modal decomposition of ^ 
in Eq. (49). In order to isolate it is easier to work with the polarization currents, which 
from Eq. (47) are given by: 


X(o})(y + x(o}W) 'cinc 


(50) 


where ^ (co) = —l/x(o}) and the inverse of J^ + x^ is guaranteed to exist (for e 7 ^ 0 , 1 ) because 
it is a Eredholm operator [172]. 

Because the extinguished power and the total LDOS are both imaginary parts of linear func¬ 
tionals of the induced fields, their modal decompositions are single sums over the resonances. 
Given the electric field from Eq. (50), the extinction of Eq. (5) and total LDOS of Eq. (15) are 
determined by overlap integrals between the incident field and the VIE-basis eigenfunctions: 



T^ext = ^ Im (E - ^ (at)) ^ 
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Because (S — (o)) ^ is a diagonal operator, it can be written as a single sum over the modes, 







simplifying the extinction and LDOS; 
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where p„ = fy • E„ Jy E„ • Einc, Pn = Ly (/y E„ • Einc.^^) , and are the VIE eigenvalues 
and the diagonal entries of E. The and p„ are normalized “oscillator strengths” representing 
the per-mode extinguished power and per-mode total density of states, respectively. Adding up 
only the oscillator strengths yields structure-independent sum rules for each quantity: 

Ep«= / lEincpdV (53a) 

n JV 

= 11 f Einc.sy • Einc.i, dv, (53b) 

fj j 

where the sum over all modes exploits the eigenbasis orthogonality, ^, and as before 

j indexes the polarization of the dipole emitter. For each quantity, the sum of the oscillator 
strengths or p„ is given by the intensity of the electric held originally incident upon the 
volume occupied by the scatterer. 

Ideally, the sum rules for the extinction and the total LDOS would lead directly to limits, 
with the numerators in Eqs. (52a,52b) bounded above by the sum rules and the denominators 
bounded below by lm^{co) = \mx{o})/\x{o})'\^, which is nonzero due to material losses. Be¬ 
cause our system is non-Hermitian, such an argument is not valid. If we had a Hermitian system 
with a conjugated orthogonality relationship [E* ■ E^ = 0, then we would have used instead 
of in Eq. (49) and the resulting amplitudes p„ and p„ would have been real and positive. 
Due to radiation losses (though not metal losses), we have complex p„ and p„ with possibly 
negative real parts, and hence it is possible to have e.g. |p„| ^ |LPn|- Such a response is a 
general feature of nonnormal dynamics, where there can be signihcant amplification beyond 
what one would expect from the resonances, as the pseudospectral level curves [140] may look 
very different from typical circles centered at the eigenvalues. 

Energy conservation prevents such responses from surpassing the sum rule limits that are 
obtained when all of the oscillator strength in Eqs. (53a,53b) are concentrated at a single res¬ 
onance (with the caveat that Einc.iy —E*|j. ^although for the dominant 1 /r^ quasistatic term 
this makes no difference). The optimal resonance is given by 

^ext.opt — Rc^ (®) (54) 

which is on the real line, as close to the material parameter ^ (w) as possible. This choice of 
eigenvalue leads to the extinction limit in Eq. (32b). Imposing energy conservation on the ab¬ 
sorption, scattering, and radiative and nonradiative LDOS integrals in the VIE approach yields 
the limits of Eqs. (32a,32b,34a,34b), but we will not prove that here. 

An interesting possibility that arises from the VIE approach is the potential existence of 
further sum rules. In addition to considering the sum of oscillator strengths in Eq. (53a), one 




can consider the sum of eigenvalues, weighted by the oscillator strengths; 


Y^^nPn 

n 

= '£[ /Er(x)E"(x)^„E"(x')Ef (xO 

jj J J 

= LJ J /Er(x)Ej'(x)S,-,(x',x")E^(x")Ef(x') 

= 1J Er(x)^;Kx',x)Ef (x') (55) 

where we used the resolution of the indentity ^„E"(x)E^(x") = 5jjt5(x —x") to simplify the 
third line. In the surface-integral representation of quasistatic electromagnetism it can be shown 
the oscillator strengths and relevant eigenvalues for extinction averaged over all angles is con¬ 
strained to satisfy Y.n^nPn/T.nPn = 1/3, reducing the possibility for all-angle response rela¬ 
tive to the single-angle response [45,47]. Equation (55) and its higher-order counterparts (e.g. 
Ln ^nPn ~ may yield stricter sum rules under various incident fields, reducing the 

possible response. 

B. Bounds on the 0(kL) term in the LDOS limits 

The bounds on the radiative and nonradiative LDOS, Eqs. (34a,34b), take into account the 
1/r^, 1/r^, and 1/r terms in the free-space dyadic Green’s function. Integrating the 1/r term 
over a half-space (or a spherical shell, or any other structure separated some distance d from 
the source) yields the 0(A:L) term in Eqs. (34a,34b), which diverges as the size L —> oo. As 
discussed in the text, this divergence is unphysical. It results from deriving the optimal current 
as proportional to the incident field, P ~ Ejnc, which is appropriate and feasible for the evanes¬ 
cent waves, but which for the 1/r term yields a physically impossible fixed energy density over 
infinite space within a lossy medium. 

Despite the divergence of this term, for finite object sizes its contribution to the limits of 
Eqs. (34a,34b) is actually very small. One does not even have to consider a finite object but a 
finite interaction distance: L represents the largest distance over which polarization currents in 
the metal generate nonzero scattered fields at the dipole source. Even wavelength-scale lengths 
are upper bounds to reasonable interaction distances in a lossy medium. Table 1 compares the 
near-held limit in Eq. (35a) to the full limit with all terms in Eq. (34a) as well as the full limit 
without the O {klfij term. 


PO “ Imj 

We take |iJ(P/Imx = 1 for simplicity. One can see from Table 1 that the near-held limit is a 
very good approximation to the overall limit for realistic separations, and that the far-held term 
contributes at most 0.03% for the cases considered. 

C. Suppressing absorption 

As discussed in the text, the optimal scattering limits are reached when absorption and scatte¬ 
ring are exactly equal. Eor some applications (e.g. solar cells enhanced by plasmonic particle 
scatterers [4]), however, parasitic absorption is detrimental and needs to be avoided, even if the 
per-volume scattering is reduced. Here we present alternate limits to Eq. (32a) to account for 
absorption suppression. 
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Table 1. Tabulation of higher-order term contributions to radiative LDOS limits. 


kd 

kL 

Eq. (35a) 

Eq. (56) 

Rel. Error 

Eq. (34a) 

Rel. Error 

0.0628 

1 

125.98 

127.97 

1.56% 

128.02 

h-0.03% 

0.01 

1 

31250 

31257.25 

0.0232% 

31257.3 

H-0.0001% 

0.01 

10 

31250 

31257.25 

0.0232% 

31257.68 

h-0.0013% 


We define a fraction / that is the ratio of absorption to extinction, 

/=—— -. (57) 

Pahs + T’scat 

Suppose we define our figure of merit as the maximum scattering cross-section per unit volume 
subject to the condition that / is smaller than some maximum ratio /max- 


max fscat 

S.t. f < /max- 


As in Sec. 3, we use standard Lagrangian optimization techniques and the the fact that the 
constraint on / is active (/ = /max) by the KKT conditions. By the same steps as in Sec. 3, it is 
straightforward to show that the per-volume scattering cross-section is limited by: 


t^scat 

~v~ 


/max(i 


(58a) 


Eq. (58a) is maximized at /max = 0.5, corresponding to the unconstrained optimum in Eq. (29a). 
Eor significantly reduced absorption (say < 2%), there is a significant penalty (> lOx) to the 
maximum per-volume scattering. 


D. Quasistatic cross-section of a coated sphere 

Also shown in Eig. 3 are the cross-sections of coated spheres that have been optimized by the 
same procedure as the ellipsoids. Coated spheres, with dielectric cores and metallic shells, are 
a common structure for tunable resonances across visible and infrared frequencies [31, 110 , 
173-175], but one can see that their response falls short of the limits. To understand why the 
performance falls short, we consider the quasistatic absorption cross-section, which is known 
analytically [32] and can be compared to Eq. (39). Eor simplicity we assume the particle core 
has the same permittivity as the shell, which has only a small effect on the response but enables 
us to write the typical [32] coated-sphere cross-section ffcs in the form 
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where fy is the metal volume fraction, Lq = 1/2—\ /6-\/1 + 8 /y and L\ = 1/2 -b 1 /()^/\ + Sfy 
are structural depolarization factors, and as in Sec. 3 we define ^ (o) = —l/x{(o)- One can see 
that there are two quasistatic resonances that arise from coupling the plasmons at the interior 
and exterior interfaces. At visible and infrared frequencies, the material parameters of typical 
metals satisfy | Rejl ^ 1 , such that a single resonance dominates the response and the volume 
fraction fy of the metal must be small (as in Eig. 3). The cross-section is then given by 
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V y /cs 3 Im^(a)) 3 lmx{o)) 


(60) 
























where we see that the absorption cross-section is proportional to |Rex|/Imx instead of 
IXp/Imj due to the small optimal volume fraction, which decreases as Rej increases. There 
is an additional 33% reduction due to the second “material” resonance [45] at Li. Additionally, 
tuning the resonant wavelengths requires very large core-shell thickness ratios, roughly 25:1 at 
X = lOOOnm, whereas at the same wavelength the optimal ellipsoids have aspect ratios of less 
than 10:1. 


E. Propagation length of a plasmon near the surface-plasmon frequency 

It is sometimes stated [33,176-179] that for any frequency at which | Ree,„| ^ Ime^, for metal¬ 
lic permittivity e„,, the complex surface-parallel wavevector for a plasmon at a metal-dielectric 
interface is given by: 


® / £'m£d ^ . £m f £'m£d 

c e'm + ^d 2{e;„)^\e!„ + £dJ 


( 61 ) 


where is the dielectric permittivity and we use single and double primes to denote real and 
imaginary parts, resp. However, this expression is only an approximation, and it is only accurate 
at low frequencies where the plasmon confinement to the surface is not highly subwavelength. 
In particular, it is not accurate near the resonant frequency at which e',« —£ci, even if |e;i»e". 
The frequency-dependent propagation length of the plasmon is given by: 

which according to Eq. (61) would go to zero at e', = —Ed, even for a lossy medium (as is 
claimed in Ref. [177]). This is not correct. 

The exact dispersion relation for a surface plasmon is [33] 


^ t Efj^Ed 
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(63) 


The approximation in Eq. (61) arises from a Taylor expansion of the imaginary part of £,„ in 
the denominator of Eq. (63). However, such a Taylor expansion is only valid if | -f Cj | ^ e". 
Thus we see that this approximation fails near the surface-plasmon frequency, which is the 
exact regime where significant subwavelength confinement is possible. If we define the surface 
plasmon frequency CUsp such that e',(a)sp) = —£d, then the dispersion relation is approximately 
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assuming |e',| ^ e" (as above). It is then straightforward to derive the propagation length at 
CO = cUsp: 
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(65) 


which does not go to zero except in the case of zero loss. (In the case of zero loss the group 
velocity also goes to zero, with the same dependence on e", such that the decay time is nonzero 
and finite. As the material loss goes to zero, the energy loss per unit time remains nonzero 
but finite, because the resonant field enhancement balances the vanishing material loss.) Thus 
near the surface-plasmon frequency the propagation length is proportional to yirnCm, not 
(Ree)^/Ime (which is the scaling of Eq. (61) and which only applies in the low-frequency 
Sommerfeld-Zenneck regime). 















F. Optimal LDOS of a thin film 

We derive Eq. (43), the optimal LDOS of a thin film at frequencies CO ^ (Usp. We measure the 
LDOS a distance d from the thin film, which we take to have a thickness t. The LDOS of a thin 
film, relative to the free-space LDOS po, is [69] 

^ ^ j (Im52i) dkp (66) 

where kp is the magnitude of the surface-parallel wavevector, and we have assumed the primary 
contribution to the integral comes from high-wavevector (p-polarized) waves for which « 
^-ikpd^ For small imaginary permittivity Ime ^ |Re£|, which is a requirement to access the 
large wavevectors, lm52i is a sharply peaked function centered around a parallel wavevector 
kpo that depends on the thin-film thickness t by [95] 

kpo « ^ (67) 

which we will use below to determine the optimal thickness. The integral can be approximated 
by 

j (Im52i) dkp « [Im52i]„ax Akp (68) 

where we have assumed that Im52i is a Lorentzian with full-width half-max of Akp. As derived 
in Ref. [95], the “reflectivity-bandwidth product” is 


~ 2kp0 


(69) 


Thus the LDOS enhancement is proportional to k^^e which reaches a maximum for a 

thickness such that kpo = 3/2z, yielding an optimal LDOS of 
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(70) 


G. 

The 


Frobenius norm of the dyadic Green’s function 

dyadic Green’s function defined in Eq. (33) is given by [51] 
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where a = kr. The Frobenius norm is given by (momentarily denoting complex con¬ 

jugation with an overline) and can be computed by use of the identities (repeated indices indi¬ 
cating summation) 


5ij5ij = 5ii = 3 

XiXj XiXj 


XiXj _ XiXi _ 

'o' ^2 “ ^ 


= 1 


which yield the expression in Eq. (33). 
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